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Depression Disorder Integrated Model (DDIM) 

A method of choosing treatment based on selection (after outcome predictions 
are made) and then a method for generating a profile for patient monitoring 
patient response to specific treatment based on a method for prediction 
treatment specific response patterns. 



Diagnosis and other baseline data are collected and transmitted as input to 
the "Outcome Predictor." The outcome predictor has been trained on other 
(prior) patients baseline and outcome data. Predictions and confidence limits 
are generated for a number of treatment options. 

Treatment selection suggestion is made by the system for evaluation and 
consideration by a clinician, who has the final responsibility for the 
treatment selection. 

The baseline data and the chosen treatment are input to the Pattern 
Predictor 
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Training Cycle for Each Treatment Group 
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Figure 3-2: Training cycle for each treatment group. Each of the two models with separate sets 
of parameters, but with a same architecture, was trained by cycling through individual data within 
the respective treatment group. After each cycle, the cost function (which reflects the degree of fit 
of the model predictions to the actual data) was evaluated to determine the completion of training. 
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Figure 3*3: Overview of the recovery model. Ellipses represent symptom factors, arrows 
between the symptom factors indicate each symptom factor can influence every other symptom 
I actor . 
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Figure 3*4: The annotated second order differential equation used to model the pattern of 
recovery. 
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Figure 3-5: Direct effects of treatment are either immediate (step function) or delayed (sig- 
moid function). Delays are estimated by treatment from the patient data, using an optimization 

procedure. 
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Figure 3-6: Direct effects and interactions in 



the recovery model. u, represents the strength 



of the immediate effect of treatment on symptom node i; „ represents the strength of the de ayed 
effect of treatment on symptom i; and and represent the interaction between the symptoms: 
^ length of the effect of symptom i on symptom," and the strength of the effect of symptom j 

on symptomi, respectively. 
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Figure 3-7: Overview of Training Process. Parameters were initialized with regression matrix 
which was ca.cu.ated from actual symptom values (ASV) by correlation and regression analyst 
The mo de used these miti-1 parameters to predict symptom values (model symptom values, MSV) 
of each pat.ent from baseline. The optim.zaUon process iteratively modified parameters to minimiz 

the discrepancy between MSV and ASV. 
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Figure 3-8: A schematic description of the cost function L. The formula inside 
the integration has two terms. By minimizing the first term, discrepancies between 
estimated and actual patterns of recovery are minimized. The second term is a 
constraint term which states that the differential equation must hold. 
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(b) CBT patient 195 
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(c) DMI patient 155 
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Figure 3-10: Individual patterns of recovery. Plots, except the one at the bottom right, 

show patterns of recovery. Solid lines show actual patterns, dotted lines show predicted patterns. 
Numbers shown at vertical axis are scaled such that the possible maximum symptom factor value 
yields 1.0. The plot at the bottom right shows the error (I) on the ordinate axis plotted against 
number of training cycles on the abscissa. 
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(a) Predicted patterns of recovery by shunting equations 
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(b) Predicted patterns of recovery by second order equations 
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(a) Predicted and actual mean patterns of recovery (CBT) 
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(b) Predicted and actual mean patterns of recovery (DMI) 

Figure 3-11: Predicted and actual recovery patterns of patient data mean, (a) mean of six 
CBT responders (b) mean of six DMI respondent In both (a) and (b), all plots except the sub- 
plots at the bottom right, show mean patterns of recovery. Solid lines show actual mean patterns, 
dotted lines show predicted mean patterns. Numbers shown at vertical axis are scaled such that 
the maximum possible symptom factor value yields 1.0. The plot at the bottom right in both (a) 
and (b) shows the error L plotted on the ordinate axis against number of training cycles on the 
abscissa. Pat 0 indicates the patient averaged patient data were used. 
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Figure 3-13: Predicted CBT and DMI temporal response sequence of symptom improvement. 
The vertical axis shows the mean half reduction time in weeks; the horizontal axis has no meaning. 
Symptom names are placed vertically at their mean half reduction time. Significant difference 
(p < 0.05 after rounding) between half reduction times of energy and early sleep disturbance. There 
is a trend (p < 0 10 after rounding) for a split between cognitions and work in CBT responders. In 
DMI there were no significant differences (or trends) in the sequence. 
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Figure 3-14: Comparison of model's predicted immediate (left) and (b) delayed (latent) direct 
effects of treatment on symptoms for Cognitive Behavioral Therapy and Desipramine. A solid line 
represents CBT coefficient values, a dashed line represents DMI coefficient values. Symptom are 
represented along the x-axis. The coefficient values the parameter optimization procedure indicate 
the strength of the effect on the symptom at the time the effect takes place, and are placed on the 
y-ax.s. For example, the delayed effect of cognitions for desipramine occurs at 3.4 weeks with a 
magnitude of almost 1.5, whereas the delayed effect of cognitions of CBT takes place at 1 2 weeks 
and has a magn.tude of about 4.2. A zero indicates that the model predicted the symptom would 
worsen initially. 
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Figure 3-15: Graphic representation of the sequence of symptom factors in recovery with 
Cognitive Behavioral Therapy treatment for the second order system. Vertical positions of the 
symptoms represent half-way-reduction time, arrows represent strong impacts and interactions, 
and corresponding numbers indicate the strength of the impact or interaction. See text for formula. 
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Figure 3-16: Graphic representation of the sequence of symptom factors in recovery with 
Desipramine treatment for the second order model. Vertical positions of the symptoms represent 
half-way-reduction time, arrows represent strong impacts and interactions, and corresponding num- 
bers indicate the strength of the impact or interaction. Dotted arrows show the interactions that 
operate in loops. See text for formula. 
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Table 4.1: Prior attempts to predict outcome from Hamilton Depression Rating Scale scores 
(symptoms) and severity (total HDRS) of depressive symptoms. Independent variables are the 
HDRS items described under the column Symptom. Numbers are indices into Table 4.2 where the 
citations are listed and indicate that the symptom at baseline was not found to be a significant pre- 
dictor of outcome. Underlined numbers are also indices into Table 4.2, but indicate the symptom at 
baseline was found to predict outcome. Blank entries indicate that the significance of the symptom 
to predict outcome was not reported. (+) indicates greater intensity predicted better response; (-) 
indicates greater intensity predicted poorer response. See Table 4.2 for key to citation references. 
Other=S-adenosyl methionine (Carney et al., 1986), ECT or unspecified antidepressant medica- 
tion (Hinrichsen & Hernandez, 1993), unspecified antidepressant medication (Katon et al., 1994); 
Tricyclic=tricyclic antidepressant medications; MAOI=monoamine oxidase inhibitor antidepres- 
sant medications; SSRI=selective serotonin reuptake inhibitors; ECT=electroconvulsive therapy; 
Ami-amitriptyline; Nor=nortriptyline; IMI=imipramine; DMI=desipramine; Clo=clomipramine; 
Lev=levoprotiline; Map=maprotiline; Phe=phenelzine; Tran=tranylcypromine; Flu=fluoxetine; 
IPT=interpersonal therapy; CBT=cognitive behavioral therapy. 
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Figure 44: Nonlinear mapping of backpropagation. Each hidden node finds a direction in the 
input space (illustrated by an arrow perpendicular to a small square piece) to which the output 
is sensitive to. The output of each hidden node goes through a nonlinear output function before 
being weighted and summed at the output node. 
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Figure 4-2: Schematic representation of the effect of normalizing transformations on reducing 
nonlinearity of score-to-output relationships (or skewedness of distributions). In the transformation, 
the area under the curve is preserved. The transformation redistributes the position of the data 
values along the x-axis in order to preserve the areas under the curve between adjacent scores values 
while redistributing these data to best approximate a normal distribution. Equal areas under the 
curve between percentiles map to equal areas under the curve in the new distribution. 



